Spectral representation of the effective dielectric constant of graded composites 
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We generalize the Bergman-Milton spectral representation, originally derived for a two-component compos- 
ite, to extract the spectral density function for the effective dielectric constant of a graded composite. This 
work has been motivated by a recent study of the optical absorption spectrum of a graded metallic film (Ap- 
plied Physics Letters, 85, 94 (2004)) in which a broad surface-plasmon absorption band has been shown to be 
responsible for enhanced nonlinear optical response as well as an attractive figure of merit. It turns out that, un- 
like in the case of homogeneous constituent components, the characteristic function of a graded composite is a 
continuous function because of the continuous variation of the dielectric function within the constituent compo- 
nents. Analytic generalization to three dimensional graded composites is discussed, and numerical calculations 
of multilayer composites are given as a simple application. 

PACS numbers: 77.22.Ej, 42.65.An, 42.79.Ry 



I. INTRODUCTION 

In graded materials J3l the physical properties may vary 
continuously in space making them distinctly different from 
homogeneous materials. Hence, composite media consisting 
of graded inclusions have attracted much interest in various 
engineering applications |2j], such as reduced residual and 
thermal barrier coatings of high temperature components in 
gas turbines, surface hardening for tribological protection, and 
graded interlayers used in multilayered microelectronic and 
optoelectronic components 010]. 

Like graded materials, thin films are of great interest in 
many practical applications and often possess different opti- 
cal properties |5| in comparison to bulk materials. Recently, 
it was found experimentally that graded thin films may have 
high relative dielectric permittivity as well as a flatter tem- 
perature characteristic of permittivity |0] than single-layer 
films H. 

The traditional theories used to deal with the homogeneous 
materials however fail to deal with composites of 

graded inclusions. To treat these composites, we have recently 
developed a first-principles approa ch l lOl Hill and a differen- 
tial effective dipole approximation 11211 . 

This work has been motivated by a recent study of the op- 
tical absorption spectrum of a graded metallic film fl3ll . In 
that work, a broad surface plasmon absorption band was ob- 
served in addition to a strong Drude absorption peak at zero 
frequency. Such a broad absorption band has been shown to 
be responsible for the enhanced nonlinear optical response as 
well as an attractive figure of merit (the degree of optical ab- 
sorption). Yuen et al. iTT^ pointed out that such an absorp- 
tion spectrum, being related to the imaginary part of the ef- 



fective dielectric constant, should equally well be reflected 
in the Bergman-Milton spectral representation of the effective 
dielectric constant fillllql . 

Bergman-Milton spectral representation was originally de- 
veloped for calculating the effective dielectric constant and 
other response functions of two-component composites fl5l 
[kJ. However, the two concerned components are all homo- 
geneous. Therefore, it is worth extending the spectral repre- 
sentation to graded composite materials. The work on graded 
films is just a simple example of a more general graded com- 
posite in three dimensions. One of the main purpose of this 
work help to identify the physical origin of the broad absorp- 
tion band. It turns out that, unlike in the case of homogeneous 
materials, the characteristic function of a graded composite is 
a continuous function because of the continuous variation of 
the dielectric function within the constituent component. 

Moreover, we apply our theory to a special case of graded 
composites, i. e., multilayer material, which is more conve- 
nient to fabricate in practice than graded material fl7ll . and 
many algorithms are now available for designing of multilayer 
coatings 1181 Il9tl . Thus, the present work is necessary in the 
sense that we shall discuss the multilayer effect as the number 
of layers inside the material increases. In this regard, this work 
should be expected to have practical relevance. As the num- 
ber of layers N increases, we shall show a gradual transition 
from sharp peaks to a broad continuous band until the graded 
composite results are recovered by the limit of N — > oo. 

The paper is organized as follows. In Sec. [II] the general 
derivation of the spectral representation for graded composites 
is presented. In Sec.|lll] we describe the model and present an- 
alytical results for the spectral representation of the effective 
dielectric constant of a graded film with an interface, as well 
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as a graded sphere. In those cases, the Bergman-Milton for- 
malism has been modified for graded composites. We further 
obtain an analytic form for the spectral density function of a 
multilayer film and a multilayer sphere in Sec. llVI Numerical 
results are presented in Sec.|Vj and discussion and conclusion 
are given in Sec. lVIl 



U. FORMALISM 

We consider a two-component composite in which graded 
inclusions of dielectric constant ei(r) are embedded in a ho- 
mogeneous host medium of dielectric constant £2- It is noted 
that the dielectric constant ei(r) is a gradation profile as a 
function of the position r. And we will restrict our discus- 
sion and calculation to the quasi-static approximation, i. e., 
dc/uj < 1, where d is the characteristic size of the inclusion, 
c is the speed of light in vacuum and uj is the frequency of 
the applied field. In the quasi-static approximation, the whole 
graded structure can be regarded as an effective homogeneous 
one with effective (overall) linear dielectric constant defined 
aslEJ 
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where Eq is the applied electric field along z direction, E and 
D are the local electric field and local displacement, respec- 
tively. 

The object of the present section is to solve the Laplace's 
equation 



V • (e(r)V</>(r)) = 



(2) 



subject to the boundary condition <po = —Eqz. The dielectric 
function e(r) varies from component to component but has a 
fixed mathematical expression for a given component. It can 
be expressed as 1151 



e(r) = e 2 



1 



-r)(r) 



(3) 



where s = [1 — e rc f/e2] _1 is the material parameter and e re f 
is some reference dielectric constant in the graded component. 
The characteristic function 77 (r) is may be written in terms of 
a real function /(r) as 



7]{t) 



1 + / (r) in inclusion, 
in host, 



which accords for the microstructure of graded composites. 
The function /(r) depends on the specific variation of the di- 
electric constant in the inclusion component. For homoge- 
neous constituent component, i.e., /(r) = 0, rj(r) = 1 in 
the inclusion component, while rj(r) = in the host medium. 
For graded systems, rj(r) can be a continuous function in the 
inclusion component because of the continuous variation of 
the dielectric function within the inclusion component. Thus, 
Eq. Q can be solved 



-E z + - s J dV'r](r')V'G(r-r')- V'0(r'), 



(4) 



where G(r — r') is a Green's function satisfying: 

V 2 G(r-r') = -<S 3 (r-r') forrinV, 
G = for r on the boundary. 

In order to obtain a solution for Eq. (0, we introduce an 
integral-differential Hermitian operator T, which satisfies 



T = J dV'r)(r')V'G(r-r')-V, 
and define an inner product as 

= / dVrj{r)X7<f)* ■ Vtp. 



(5) 



With the above definitions, Eq. (0J can be simplified to 



0(r) = -E z+-f<t>(r). 

s 

Let s n and \<f> n ) be the nth eigenvalue and eigenfunction 
of operator T. Then, the generalized eigenvalue problem be- 
comes 

V • (7,(r)V0 n ) = s„V 2 </V 
The potential \<fr) can be expanded in series of eigenfunctions, 
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where we choose Eq = 1 for convenience. Since rj(r) is 
a real function, the eigenvalues s n will be real. Also, for 
graded component, 7y(r) is a continuous function, which will 
cover the full region, i.e., —00 < rj(r) < 00. Therefore, 
the eigenvalues s n , which depend on the continuously graded 
microstructure rj(r), do not lie within the interval [0, 1] but ex- 
tend to — 00 < s n < 00 as first pointed by Gu and Gong 1I2U1 
for three-component composites case. However, eigenvalues 
s„ still lie in [0, 1] for < r?(r) < 1. 

We are now in the position to find an analytical repre- 
sentation for the effective dielectric constant e e according to 
Eq. Q. We take advantage of Green's theorem, the boundary 
condition </>o = —z, and the Maxwell equation V • D = to 
obtain the effective dielectric constant 

~ = \A (-W).DdV 

£2 £2V J 
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If we now introduce the reduced response fl^l 

F{s) = l- e f, 
and substitute Eq. (|6jl into Eq. @ we find 



(8) 
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We can now express the effective dielectric constant as 

In \ 



<~2 



where /„ is given by 



(9) 



V {(j>n\<j>n) 

Using the above equations, we obtain the following sum rule 

n 

= V J rf ^( r ) Vz ' Vz 

= iy dV^r). (10) 

It is worth noting that the sum rule will not equal to the vol- 
ume fraction of inclusion. This is different from the Bergman- 
Milton spectral representation for two homogeneous systems, 
in which the sum rule equals to the volume fraction of the 
inclusion. 





FIG. 1: (a) Dielectric profile of various graded films at h — 0.6. (b) 
Dielectric proifle of various graded sphere with unit radius. Parame- 
ters used: ei = 1 and s — — 2. 



When the operator f has a continuous spectrum, Eq. l|9} 
should be replaced with the integral form 



e 2 1 



ds 



, m(s') 



(11) 



where m(s') is the spectral density function. Then, the re- 
duced response becomes 



F(s) 



ds 



,m(s') 



(12) 



If we write s as s + i0 + , the right side of Eq. J 1 2b becomes 

„ f , .mis') , „ 

P I ds'— - mm(s), 



and thus, m(s') is given through the limiting process 



m(s') = Im[F(s' + iO+)]. 

7T 



(13) 



This final result is identical in form to Bergman's expres- 
sion for the analogous function in scalar composite materials. 
However, there are differences in the derivation, namely, the 
definition of the inner product Eq. l|5}, the continuous graded 
microstructure r?(r), the sum rule, as well as the range of 
eigenvalues s n . 

From Eq. (lilt it is evident that if the spectral density func- 
tion m(s') is known, the effective dielectric constant can be 
obtained accurately, and vice versa. The spectral represen- 
tation has been used to analyze the effective dielect ric p rop- 
erties of composites. Recently, Levy and Bergman 12211 also 
used it in their study of nonlinear optical susceptibility. In this 
regard, Sheng and coworkers l23ll developed a practical algo- 
rithm for calculating the effective dielectric constants based 
on the spectral representation. In what follows, we restrict 
ourselves to a graded composite both in one dimension and 
three dimensions, as well as corresponding multilayer com- 
posites. 



in. SPECTRAL DENSITY FUNCTION OF GRADED 
COMPOSITES 

A. Spectral density function of a graded film 

We consider a graded dielectric film of width L , in which 
two media meet at a planar interface as shown in Fig. ^( a )- 
The first medium e\(z) varies along z— axis, while the second 
medium e 2 is homogeneous. We define the graded microstruc- 
ture as 



r)(z) = 



1 + az < z < h, 
h < z < L, 



(14) 



where a and h are real constants. They can be varied to de- 
scribe different graded films. Thus, according to Eq. 0, the 
dielectric function of graded film can be expressed as 



e{z) = e a 1 - 



(15) 
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Owing to the simple geometry of a graded film, we can use 
the equivalent capacitance of a series combination to calculate 
the effective dielectric constant as 



e(z) 



dz. 



(16) 



Substituting Eqs. (I14> and Jl 51 into Eq. d!6i . we obtain 




FIG. 2: (a) Spectral density function of a graded film without an 
interface, i.e., h = 1.0. (b) Spectral density function of a graded 
film meeting a homogeneous medium at an interface h = 0.5, and 
e 2 = 1. 
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with the assumption L = 1, 

We are now in a position to extend the Bergman-Milton 
spectral representation of the effective dielectric constant fl5l 
Hal to a graded film. For a graded system, 77(2) can 
be a continuous function in the inclusion medium. Using 
Eqs. il 1> — ( fT3t . we obtain the spectral density function for 
a graded film as 



as arg 



m(s') 



s-1 



s — ah — 1 



s'arg 



s-1 



s — ah — 1 



+ \a(h- 1) - s'ln 



s'-l 



s/ — aft, — 1 



where s' = Re[s] and arg[- • • ] denote the arguments of com- crostructure as 
plex functions. 



B. Spectral density function of a graded sphere 

The above theory can be generalized to graded composites 
in three dimensions. We consider a graded sphere with di- 
electric constant ei(r) embedded into a homogeneous host 
medium with dielectric constant €2- The dielectric constant 
of the graded sphere ei(r) varies along the radius r. We can 
obtain the effective dielectric constant of a graded sphere us- 
ing the spectral representation. We consider the graded mi- 



rj(r) 



1 + ar < r < R, 







r > R, 



where R is the radius of the graded sphere. Thus, from Eq. 
the dielectric constant in the graded sphere is given by 



<r) = e 2 1 



(17) 



In the dilute limit the effective dielectric constant of a small 
volume fraction p of grad ed spheres embedded in a host 
medium is given by l24tl25ll 



e e = e 2 + 3e 2 pb, 



(18) 
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where b is the dipole factor of graded spheres embedded in 
a host as given in Ref. fioll . Using Eq. and Eq. dl St . the 
reduced response can be obtained as 



F(s) = ~3e 2 pb. 



(19) 



Thus, the spectral density function of a graded sphere can be 
given through a numerical evaluation of Eq. Jl 3I >. 



IV. SPECTRAL DENSITY FUNCTION OF MULTILAYER 
COMPOSITES 

A multilayer composite is a special case of graded com- 
posites. The gradation becomes continuous as the number of 
layers approaches infinity. To investigate the multilayer effect, 
we shall use a finite difference approximation for the graded 
profile (Eqs. Jl 5i and (I17» for a finite number of layers. To 
mimic a multilayer system, we divide the interval [0, 1] into N 
equally spaced sub-intervals, [0, z{\, (zi, Z2), (zn — 1, 1]. 
Then we adopt the midpoint value of e(z) for each sub-interval 
as the dielectric constant of that sublayer. In this way, we cal- 
culate the effective dielectric constant, eigenvalues, as well 
spectral density function for each TV. It is worth noting that 
the results of TV — > 00 (e. g., TV = 1024 recovers the results 
of graded composites. 

In addition to multilayer films, we can use the above ap- 
proach to study the much simpler problem of a two-layer film. 
In this system, we have two layers of dielectric constants e\, 
e 2 , and host e . Thickness are hy, h(l — y), and 1 — h, respec- 
tively, where y is the length ratio between component ei and 
component e 2 . We also define two microstructure parameters, 
771 and 772. If we let s = 1/(1 — ei/eo), then 771 = 1, and 
V2 = (eo — e 2)/(eo — According to Eq. il5i , the effective 
dielectric constant of the two-layer film is now given by 



1 



hy + M 1 - v) 

ei e 2 



1 



e e ei e 2 e 

According to Eq. the reduced response can be given by 



F(s) 



F] 



s — Si s — S2 



(20) 



where 

Fx -- 

F 2 = 
si 

S2 



K s i{y - vv + v) - y) 

si - s 2 

H s 2(y - yv + v) - v) 

Si - s 2 

= o I 1 - h (y ~ yy + y) + y 



-Vi^pl + h) + {l- h(y -yr] + i]) + rif 

\ [1 - Kv - yy + y) + v 



+ x /4t7(— 1 + h) + {l- h(y - yr] + rj) + r)) 2 



sum rule expressed by Eq. dlOi is obeyed. It should also be 
noted that there are two poles in the expression for the reduced 
response corresponding to two peaks in the spectral density 
function. If h = 1, then si = 0, that is, one peak is located at 
zero, which is explicitly shown in Fig. [2 a). 




FIG. 3: (Spectral density function of a graded sphere with volume 
fraction p — 0.1. 

Similarly, we can also apply our graded spectral represen- 
tation to a single-shell sphere of core dielectric constant ei, 
covered by a shell of e 2 , and suspended in a host of eo. In 
this example, we can also define two microstructure parame- 
ters r/i and rj 2 . If we let s = 1/(1 — ei/eo), then ?/i = 1, and 
y-2 = (eo — £2) / (en — ei). The dipole factor of single-shell 
sphere is given I24l 12511 



b = 



£2 - e + (eg + 2e 2 )xf 3 
e 2 + 2e + 2(e 2 - e )xf 3 ' 



where / is the ratio between radius core and radius shell, and 
x is given by 

ei - £2 



ei + 2e 2 



Then, we can also write Eq. dl9l similarly to Eq. i20l . where 
the residues and eigenvalues are given by, 



Fi 
F 2 

S2 



-3psi[(-l + r))y 3 - 77] - r)p[l - 2(-l + r/)y 3 + 2t]} 
3(si - s 2 ) 

3ps 2 [(-l + r))y 3 - rj\ + r)p[l - 2(-l + i])y 3 + 2rj\ 



3(si - s 2 ) 



1 + 3?7 - y/i + (2- 82/ 3 )?7 + (1 + 8y 3 )77 2 



1 + 3?7 + y/i + (2- 8ty 3 )r/ + (1 + 8y 3 )r 



From the sums of F\ and F 2 and the integral of graded mi- 
crostructure rj(z) given by Eq. dl4> . we can check that the 



Analysis shows that the spectral representation for TV = 2 
contains two simple poles corresponding to two peaks in the 
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FIG. 4: Spectral density function of various multilayer film with eo = 1 and a = —0.8. 



spectral density function. Therefore, we draw the conclusion 
that, N peaks are a result of N layers. Moreover, N — 1 
peaks will accumulate into a continuous broad absorption 
spectrum when N tends toward infinity, which can be seen 
from Fig. Elf) and0f). 

V. NUMERICAL RESULTS 

We are now in a position to do some numerical calculations 
of the spectral density function from Eqs. (lilt and (1131 . A 
small but finite imaginary part in the complex parameter has 
been used in the calculations. Without any loss of generality, 
we choose L = 1 and R = 1 for convenience. We show 
the effect of different graded profiles, as well as the effect of 
the thickness of the inclusion. It should be noted, that in all 
figures the range of s is limited to [0,1], because we chose 
— 1 < a < which limits the value of r\ into [0, 1]. 

Fig. ^ displays the dielectric profile of a graded film 
(Fig. E a )) and a graded sphere (Fig. [flb)). This figure ob- 
viously shows that the dielectric constant varies with the po- 
sition in inclusion while a constant in host medium. Also, 
different values of a accord with different graded materials. 

In Fig.[5Ja), we plot the spectral density function m(s) of a 
graded film without an interface against the spectral parame- 
ter for various graded microstructures r){z). It is evident that 
there is always a broad continuous band in the spectral density 
function. Both the strength as well as the width of the contin- 
uous part of m(s) increase with the gradient of the dielec- 
tric profile. Thus, the previous results of the broad surface- 



plasmon band can be expected. Note that there is a sharp peak 
at s = 0, which is also present in a homogeneous film. In 
Fig-Gib), we plot the spectral density function of a graded film 
meeting a homogeneous medium at an interface for various 
graded microstructure r/(z). Again, there is always a broad 
continuous band in the spectral density function. However, 
the sharp peak has now shifted to a finite value of s, which is 
also present in a homogeneous film. 



In Fig. |3 the spectral density function of graded sphere is 
displayed for a volume fraction p = 0.1. In this case, the inter- 
face always exists. It is clear that a broad continuous function 
in the spectral density function is always observed, as well as 
the shift of the sharp peak. However, the decrease of the broad 
continuous function is more abrupt for graded sphere than for 
graded film with increasing s. 



Figures E] and [5] display the spectral density function for 
a multilayer film and a sphere, respectively. It is clear that 
there are always N sharp peaks for N layers. Moreover, it is 
worth noting that there occurs a transition from sharp peaks 
to a broad continuous band with increasing N (see Fig. Elf) 
and Fig. |3 f)), that is, the graded results are recovered by the 
limit results of N — > oo. In particular, we had obtained the 
analytical expression of spectral density function for N = 2. 
There are two resonances corresponding to the two peaks in 
Fig.Ela)andFig.|3a). 
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FIG. 5: Spectral density function of various multi-shell sphere with eo = 1, a = —0.8 andp = 0.1. 



VI. DISCUSSION AND CONCLUSION 

We have investigated a graded composite film and a sphere 
by means of the Bergman-Milton spectral representation. It 
has been shown that the spectral density function can be ob- 
tained analytically for a graded system. However, unlike in 
the case of homogeneous constituent components, the charac- 
teristic function is a continuous function due to the presence 
of gradation. Moreover, the derivation as well as some salient 
properties, namely, the sum rule, the definition of inner prod- 
uct, the definition of the integral-differential operators, and the 
range of spectral parameters, do change because of the contin- 
uous variation of the dielectric profile within the constituent 
components. It should be noted that in graded composite, the 
eigenvalues are not limited to [0, 1], and they can be extended 
to — oo < s n < oo for the full region rj, i. e., — oo < 77 < 00. 
In this work, however for simplicity, we investigated the spec- 
tral density function in < s < 1 by choosing — 1 < a < 
to limit the value of ?/ into [0, 1]. 

We also study multilayer composites and calculated the 
spectral density function versus the number of layers, to ex- 



plicitly demonstrate that the broad continuous spectrum arises 
from the accumulation of poles when the number of lay- 
ers tends to infinity. This finding coincides with the broad 
surface-plasmon absorption band associated with the optical 
properties of graded composites. 

To sum up, we have investigated the spectral density func- 
tion of graded film and graded sphere, as well as multilayer 
cases. There is always a broad continuous function in the 
spectral density function in graded composite, but simple 
poles in multilayer composite, the number of pole depends on 
the number of layers. Moreover, there is a gradual transition 
from sharp peaks to a broad continuous band until the graded 
composite results recover in the limit of N — > 00. 
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